Uniqueness of the solution to the 
Vlasov-Poisson system with bounded density 



Gregoire Loeper 
February 1, 2008 



Abstract 

In this note, we show uniqueness of weak solutions to the Vlasov- 
Poisson system on the only condition that the macroscopic density p 
defined by p(t,x) = J Rd f(t,x,£)d£ is bounded in L°°. Our proof is 
based on optimal transportation. 



1 Introduction 

The Vlasov-Poisson system (herafter (VP)) describes the evolution of cloud 
of electrons or gravitational matter through the equations 

^/+e-v S6 /+v*.v e / = o, (i) 

= ep, (2) 

where p(t,x) = f f(t,x,£)d£, and e > in the electrostatic (repulsive) case, 
e < in the gravitational (attractive) case. Here f(t,x,£) > 0, denotes the 
density of electrons (or matter) at time t e M + , position x G K 3 , velocity 
(GK 3 . Equation means 

*(t,x) = -e/ p(t,y)—^ rdy. (3) 

J R3 4vr|x - y\ 

We denote by .M(M 6 ) (resp. A / i + (M 6 )) the set of bounded (resp. bounded 
and positive) measures on R 6 . Given an initial datum /° e A / i + (M 6 ), we look 
for solutions to (JU 12)) such that 

f\t=o = f- (4) 
For T > 0, we will call / a solution to (JU ij) in £>'([0, T[xl 6 ), if 
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- feC([0,T[,M + (R 6 )-w*), 

- V<p G C c °°([0,T[xi? 6 ), 

/ f(dt<p + £ ■ V x <p + V** • V^) = - / /V|t=o, (5) 
J Jr 6 

- * solves ©. 

We will not discuss the conditions needed on / to give sense to the product 
/V-j^ or to the singular integral Q, since we will only consider the case 
where p G L°°. In this case, will be continuous, and the product fV x ^ 
will be well defined for / a bounded measure. 
Our result is the following: 

Theorem 1.1. Given f° in M + (R 6 ), given T > 0, there exists at most one 
weak solution to 0) in V'([0, T[xR 6 ) such that 

sup HpHloo^j) < +oo. (6) 
te[o,T[ 

Remark 1. To establish the existence of a solution to (VP) satisfying 
the bound (JHJ) requires much more assumptions on the initial datum than 
what we need here ! This question is treated in 0j. 

Remark 2. Note that we do not ask for any bound on the moments of 
/, and also that we do not ask the energy to be finite. 

A sufficient condition for uniqueness had been given by Lions and Perthame 
in [3], relying on Lipschitz bounds on the initial data f°, but they expected 
a uniqueness result under the weaker assumption of bounded density. The 
Lipschitz condition had indeed later been relaxed by Robert in jH] down to 
/ G L°° compactly supported in x and £ for t G [0,T]. Here we relax the 
bound on the support of /, and we do not ask either / to be bounded in 
L°°. We only need a L°°([0, T[xIR 3 ) bound on p(t,x). Hence our result 
applies also to monokinetic solutions of ((TJ |2j). In that that case, we have 
f(t, x,£) = p(t, x)8(£ — v(t, x)) for some vector field v, and this gives formally 
a solution to the Euler-Poisson system 

d t p + V ■ (H = 0, (7) 
d t (pv) + V ■ (pv ® v) = pW, (8) 
A\& = ep. (9) 



Our proof will rely on optimal transportation, and the next section is 
devoted to recall some facts concerning this subject. A complete reference on 
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this topic is jU] . The technique we will use adapts to many similar problems, 
where a transport equation and an elliptic equation are coupled. The velocity 
field is the gradient of a potential satisfying an elliptic equation whose right 
hand side depends smoothly on the density. This has been observed in jH] in 
the case of the 2-d incompressible Euler equations and the semi-geostrophic 
equations. It is interesting to notice that our technique gives a new proof 
of the uniqueness part in Youdovich's Theorem, while the technique used by 
Robert in [H] was and adaptation of Youdovich's original proof (see [10]). 

2 Preliminary results on optimal transporta- 
tion and Wasserstein distances 

Definition 2.1. Let p\,pi he two probability measures on M. d . We define the 
Wasserstein distance between p\ and p 2 , that we denote W^Cpi, P2), by 



where the infimum runs over probability measures it on IB^xR^ with marginals 
Pi and p2- 

Then we gather several results of optimal tranpsortations in the following 
theorem. These results can be found in Benamou & Brenier j2] , McCann jHj , 
Gangbo & McCann 

Theorem 2.2 (Benamou, Brenier, Gangbo, McCann). Let p\, p 2 be 

two probability measures on M. d , such that p\,pi are absolutely continuous 
with respect to the Lebesgue measure. Assume that W2(pi,pz) < +oo. 

1. There exists a pair (p e ,v e ),8 e [1,2] such that p e E C([l, 2]; M (R d ) - 
w*) and vg is a dps measurable vector field, that satisfies 



2. For this 'optimal' path we have also, when both pi and p<i have densities 
in L°° with respect to the Lebesgue measure, 

V6» G [1,2], | l/ty 1 1 2,°° < max{||pi|| L oc, ||p 2 |U°°}- 

Remark. The path {pe,9 e [1,2]} is the geodesic linking p 1 to p 2 , with 
respect to the Wasserstein metric (see jH] , and also |Zj where this notion was 
introduced). 




d e pe + V ■ (p e v e ) = 0, 
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3 Proof of Theorem 11.11 

From now on, we assume for simplicity that j R6 f°(x,£) = 1, e = 1, and the 
reader can check that this choice does not play any role in the proof. In 
particular, the result of the previous section adapt with minor changes to the 
case of two positive measures of equal total mass. 

3.1 H 1 estimates along geodesies 

In this section, we show that the H^ 1 (M 3 ) norm of the difference p% — p 2 can 
be controlled by the Wasserstein distance between p\ and p 2 , provided both 
measures have densities in L°° with respect to the Lebesgue measure. This 
is the crucial estimate at the core of our result. 

Proposition 3.1. Let p\, p 2 be two probability measures on lR d with L°° den- 
sities with respect to the Lebesgue measure . Let ^i,i = 1,2 solve 

A^i = Pl 

in the sense of Then 

||Wi - W 2 ||L2(]Rd) < [max{||pi|| L oo, ||p 2 ||L°°}] 2 W 2 (pi,P2), 

where W 2 {pi,P2) is the Wasserstein distance between p\ and p 2 given in Def- 
inition \2.1\ 

Proof of Proposition 13.11 We suppose that W 2 (pi,p 2 ) < +00 oth- 
erwise there is nothing to prove. For 9 £ [1,2], we take pg that interpolates 
between p\ and p 2 as in Theorem 12 .2\ and vg the corresponding velocity field. 
If we consider, for every 9 £ [1, 2], ^g solution of 

= p e , (10) 

then tyg interpolates between ^I^ and ^ 2 . If we differentiate ()10|) with respect 
to 9, we obtain 

Adg^g = dgpg 

= -V ■ (pgvg). 

Note that since p±,p 2 are bounded, so will be pg from Theorem 12 .2\ hence 
using that W 2 (pi,p 2 ) = J pg\vg\ 2 < +00, we have pgVg bounded in L 2 . This 
implies dg^g £ W 1 ' 2 . 
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We integrate the above equation against dg^Q to obtain 
J |V<9^ e | 2 = -J p e v e ■ Vdg^e, 

and this yields 



< 



1 1 1/2 

Pe\\ Lao{Rd) W 2 (pi,P2)- 



Using then the second point of Theorem 12.21 and integrating over 9 G [1, 2], 
this proves Proposition 13. II 

□ 



3.2 Lagrangian formulation of the Vlasov-Poisson sys- 
tem 

Given a solution of (VP) with bounded density p on [0, T[ , we consider for 
t G [0,T[ the characteristics of equation ((H), that solve the ODE 

* = (11) 
£ = V*{t,X). (12) 

Since we assume an L°°([0, T[xIR 3 ) bound on the density p, the field V 1 ^ 
classically satisfies a log-Lispchitz condition: 

Vt G [0, T[, V(ac, n) 6 t 3 x IR 3 , |ac - y| < i 

|W(t,x) - V#(t,y)| < C|x-y[log 1 1 , , (13) 

k — 2/1 

where C depends on ||p||l°o([o,t[xR 3 )- This condition is enough to define a 
Holder continuous flow 

E(t,x,v) = (X,£)(t,x,v) 

for the ODE (|11J12J1 . where (X,£) is the pair (velocity, position) at time t of 
the trajectory having (velocity, position) equal to (x, v) at time 0. 
Then we use the following Theorem, proved in P: 

Theorem 3.2. Let u(t,x) be a vector field on M. d . Consider the ODE 

7(f)=«(t,7(*)), 
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and the PDE 

d t /j,(t, x) + V ■ (fi(t, x)u(t, x)) = 0. 
Let B C M d be a Borel set. The following are equivalent: 

(a) For all x in B, there exists a unique solution to the ODE starting at x. 

(b ) Non negative measure-valued solutions to the PDE with initial data u° 
concentrated in B are unique. 

From this result, we deduce the following corollary: 

Corollary 3.3. The potential *f> being held fixed, and satisfying A\l/ G L°°([0,t]x 
M 3 ), for any f° G A4 + (R 6 ) there exists a unique weak solution to (QJ) (i.e. in 
the sense of $)) with initial datum f° which is given by 

f(t) = 5(f, 0#/°, (14) 
where H = (X,£) solves Ml\ Note also that we will have 

p(t)=X(t,;-) # f°. (15) 

We remind the reader that the measure f(t) = H(i, -, -)#/° is defined by 
f(t)(B) = /^S^ 1 (*)(£)) for all Borel subsets B of R 6 . 

Remark. This corollary does not solve the uiqueness problem, but says 
only that if we suppress the coupling between \1/ and p, there is a unique weak 
measure- valued solution to the transport equation ((TJ) , that we can represent 
with the help of characteristics. 

3.3 Final estimate 

Given an initial distribution f°(x,£) G A4 + (R 6 ) with f R6 f° — 1> we take 
two solutions (/i,/^) to (VP) with bounded density and initial datum f°. 
We have A\& = pi,i — 1,2 in the sense of (JHJ). We then consider the asso- 
ciated characteristics Si and H 2 , where for i — 1, 2, Hj = (JQ, x, £) 
and solve (fTTl [T2J) with force field V^. Note that we will have 

fiif) = Ei(t)#f°,i = 1,2. We then consider 

Q(*) = ^ / f°(x,c)\^i(t,x,0-^{t,x,0\ 2 - 
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Remark. Notice that (Si(i), S 2 (t))#/° is a probability measure on M 3 x 
M 3 , with marginals /i(t) and / 2 (t), hence by Definition 12. 11 

Q{t)>\wi{h{t),h{t)). 

(We will repeat this argument in Lemma 13.41 ) This implies in particular 
that (Q = 0) -<=>- (/i = / 2 ). Our proof will rely on an estimate on the 
Wasserstein distance between f\ and / 2 , while the proof of [Hj was obtained 
by estimating the H -1 norm of f\ — / 2 . 

Of course Q(0) = 0, and 
tM) = [ /°(x,0(2i(t,x,0-5 2 (t,x,0)-^(S 1 (t,x,0-S 2 (t,x,0) 

Clt /id>6 



[(Xi(t, x, - X 2 (t, x, 0) • a, - x, 0)] 

(V* x (t, *i(t, x, 0) - W 2 (t, X 2 (t, x, 0)) 

The second line is bounded by Q(t), and using Cauchy-Schwartz inequality, 
the third line is bounded by 

/ r \ V2 

Q 1/2 W / /°(x,oivvi/ 1 (t,x 1 (t,x,0)-vvi> 2 (t,x 2 (t,x,0)l 2 
< Q 1/2 (t)(f /°(x,e)|v* 2 ( t ,x 1 (t,x,0)-v* 2 (t,x 2 (t,x,0)l' 

+ / /°(x,0|V* 2 (t,X 1 ( t ,x,0)-V* 1 ( t) X 1 ( t ,x,0)| 

</R 6 

Hence we have 

|Q(t) < Q(t) + Q 1/2 (t) (T^t) + T 2 (t)) 1/2 , (16) 
and we will now estimate T 2 and then T\. 

For T 2 we have, using ()15|) and Proposition 13. 11 

T 2 (t) = / p 1 {t,x)\V* 1 (t,x)-VV 2 (t,x)\ 2 
< CW*( Pl (t),p 2 (t)), 



1/2 

2 " 
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where C depends on the L°° norms of pi, p 2 . Moreover, from the very defini- 
tion of the Wasserstein distance given in Definition 12.11 we have the elemen- 
tary lemma: 

Lemma 3.4. Given Xi, X 2 , pi, p 2 as above, we have 

W 2 ( Pl (t),p 2 (t))< ( [ f (x,Z)\X 1 (t,x,0-X 2 (t,x,0(' 



Proof. The proof follows immediately from Definition 12.11 by notic- 
ing that 7T = (Xi(t), X 2 (t))#f° is a probability measure on R 3 x M 3 with 
marginals pi(t) and p 2 {t). 

□ 

Hence 

W*( Pl ,p 2 ) < [ f°(x,0\X 2 (t,x,0-X 1 (t,x,0\ 2 

< / /°(x,0|2i(t,x,0-H 2 (t,x,0| 2 = 2g(t), 

and we conclude that T 2 {t) < CQ{t) for some C depending on the L°° bounds 
on pi,p 2 . 

Now, we evaluate 7\ by standard arguments, using the log-Lipschitz 
regularity of V\l/ 2 : note first that since p iy i = 1,2 are bounded in L°°, 
V^j, i = 1,2 are also bounded in L°°, hence for any C > 0, we can take T 
small enough such that ||Hj — S 2 || iO o([ 0)T ] xR 6) < C. Thus we have, for some 
other C depending on ||pi||L°°,« = 1,2, and as long as ||E;l — S 2 || i o C < |, 

Ti = [ / (x,0|VM/ 2 (t,X 1 (t,a;,0)-VvI/ 2 (t,X 2 (t,x,0)| 2 

< C[ f°(x,0(\X 1 -X 2 \ 2 log 2 1 , W,0 
Jm 6 V l A i ~~ A 2| / 

= t/ r^od^i-^piog^iXi-x.i 2 ))^^^). 

Then we use that x i— > a; log 2 x is concave for < a; < 1/e, and we can 
assume (taking T small enough) that ||Hj — S 2 ||i,oo(r 0)T | xR 6) < 1/e, therefore 
by Jensen's inequality we have 



rx(t) 
c 



/ r^oiXx-x,! 2 ^^^) log 2 / / (x,o|X!-x 2 | 2 (t,x,o 

</IR 6 J UR 6 



< 

4 

< ^g(t)io g 2 (2g(t)). 
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Combining all these bounds in (|16j) . we obtain that 

d -Q(t)<CQ(t)(l + \og- ' 



dt^ y/ ~ Q{ty 

and we conclude by standard arguments that if Q(0) = 0, Q = on [0, T[. 
This achieves the proof of Theorem 11.11 

□ 
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